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ABSTRACT 

The discovery of multiple stellar populations in globular clusters has implications for all the 
aspects of the study of these stellar systems. In this paper, by means of A'-body simulations, 
we study the evolution of binary stars in multiple-population clusters and explore the im¬ 
plications of the initial differences in the spatial distribution of different stellar populations 
for the evolution and survival of their binary stars. Our simulations show that initial differ¬ 
ences between the spatial distribution of first-generation (FG) and second-generation (SG) 
stars can leave a fingerprint in the current properties of the binary population. SG binaries 
are disrupted more efficiently than those of the FG population resulting in a global SG binary 
fraction smaller than that of the FG. As for surviving binaries, dynamical evolution produces 
a difference between the SG and the FG binary binding energy distribution with the SG pop¬ 
ulation characterized by a larger fraction of high binding energy (more bound) binaries. We 
have also studied the dependence of the binary properties on the distance from the cluster 
centre. Although the global binary fraction decreases more rapidly for the SG population, the 
local binary fraction measured in the cluster inner regions may still be dominated by SG bi¬ 
naries. The extent of the differences between the surviving FG and SG binary binding energy 
distribution also varies radially within the cluster and is larger in the cluster inner regions. 

Key words: globular clustersigeneral, stars:chemically peculiar. 


1 INTRODUCTION 

Numerous observational studies have revealed that globular clus¬ 
ters host multiple stellar populations characterized by different 
chemical (see e.g. Carretta et al. 2009a, 2009b, Gratton et al. 2012 
and references therein) and photometric properties (see e.g. Lee et 
al. 1999. Bedin et al. 2004, Ferraro et al. 2004, Piotto et al. 2007, 
2015; Milone et al. 2008, 2010, 2012; Bellini et al. 2013). The im¬ 
plications of the presence of multiple populations in globular clus¬ 
ters and the questions raised by this discovery span all the aspects 
of the astrophysics of these stellar systems including their forma¬ 
tion, chemical and dynamical evolution along with the possible role 
played by these objects in the assembly of the Galactic halo. 

A complete scenario for the formation and evolution of 
multiple-population globular clusters requires the identification of 
the source(s) of gas with chemical properties consistent with the 
observed abundance patterns and out of which second-generation 
(SG) stars might have formed along with models following the dy¬ 
namics of this gas, the SG formation history and the subsequent 
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early and long-term dynamical evolution of the multiple-population 
cluster. Different sources of gas for the SG formation have been 
suggested in the literature (see e.g. Ventura et al. 2001, Prantzos & 
Charbonnel 2006, Decressin et al. 2007, de Mink et al. 2009, Bas- 
tian et al. 2013) and a number of studies have addressed some of 
the issues concerning the origin of the observed abundance patterns 
(see e.g. Decressin et al. 2007, D’Ercole et al. 2008, 2010, 2012), 
the formation and dynamical evolution of multiple-population clus¬ 
ters (see e.g. D’Ercole et al. 2008, Bekki 2011, Vesperini et al. 
2010, 2013, Decressin et al. 2008, 2010). The theoretical study of 
multiple-population clusters is still in its infancy and most ques¬ 
tions and challenges are just now beginning to be addressed. 

In D’Ercole et al. (2008), we have modelled the formation and 
dynamical evolution of the SG population forming from the ejecta 
of first-generation (FG) intermediate-mass AGB stars; one of the 
predictions of our study was that the SG population would form 
in a compact subsystem concentrated in the inner regions of the 
FG cluster. In a subsequent study on the long-term evolution of 
multiple-population clusters (Vesperini et al. 2013) we have shown 
that some memory of the initial spatial segregation of the SG pop¬ 
ulation should be preserved in many Galactic clusters and, indeed, 
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several observational studies (see e.g. Sollima et al. 2007, Bellini et 
al. 2009, Lardo et al. 2011, Milone et al. 2012, Beccari et al. 2013, 
Johnson & Pilachowski 2012, Cordero et al. 2014, Kucinskas et al. 
2014) have found clusters in which SG stars are more spatially con¬ 
centrated than FG stars (see also Bellazzini et al. 2012 and Richer 
et al. 2013 for interesting studies on the kinematical properties of 
FG and SG stars). 

One of the dynamical implications of the SG spatial segre¬ 
gation concerns the evolution of binary stars. Until SG and FG 
stars are completely mixed, SG stars will evolve in a denser envi¬ 
ronment where disruption and evolution of SG binaries will occur 
more rapidly than for FG binaries. In Vesperini et al. (2011) we 
have carried out an initial investigation of this issue by combining 
iV-body and semi-analytical calculations of the disruption of bi¬ 
naries in multiple-population clusters. Our study showed that SG 
binary disruption rate can indeed be significantly larger than that of 
FG binaries and that the overall binary disruption rate in multiple- 
population clusters is larger than that of a single population cluster 
without the central SG dense subsystem. 

An early investigation by D’Orazi et al. (2010), although 
based on a small sample of Ba stars, suggests that the fraction of 
SG binaries might indeed be smaller than that of FG binaries. 

As for the connection with observational studies, as already 
pointed out in Vesperini et al. (2011) and further discussed in more 
detail in this paper, it is important to emphasize the differences be¬ 
tween the theoretical predictions concerning the global binary frac¬ 
tion and the local binary fraction measured at different distances 
from the cluster centre. Specifically, we will show in this paper that 
even if the global SG binary fraction is smaller than the FG one, the 
local SG binary fraction measured in the cluster innermost regions 
may still be larger than the FG one. 

In this paper we extend the initial investigation by Vesperini 
et al. (2011) by carrying out a survey of iV-body simulations and 
including in a self-consistent way a population of binary stars in¬ 
stead of just relying on analytical estimates for binary disruption 
rates due to single/binary interactions. The iV-body simulations 
presented here, while still idealized in several aspects, allow us to 
follow in much more detail the evolution of the cluster binary pop¬ 
ulation and its dynamics. The questions we address in this paper 
concern the radial variation of the fraction of FG and SG binary 
stars resulting from the combined effect of disruption and tendency 
of binary stars to segregate in the cluster inner regions, the evolu¬ 
tion of the global binary fraction and its relation with the local frac¬ 
tion measured at different distances from the cluster centre, and the 
evolution of the binding energy of the surviving binaries. 

The results of our investigation show that differences in the 
initial structural properties of the FG and SG populations can leave 
a significant fingerprint in the properties of FG and SG binaries. 

The outline of this paper is the following. In Section 2 we 
describe the method and the initial conditions adopted in our sim¬ 
ulations. In section 3 we present our results, and in section 4 we 
summarize the main conclusions of our investigation. 


2 METHODS AND INITIAL CONDITIONS 

The iV-body simulations presented in this paper have been car¬ 
ried out with the GPU-accelerated version of the code NBODY6 
(Aarseth 2003, Nitadori & Aarseth 2012) and run on the BIG RED 
II cluster at Indiana University. 

In the initial conditions we have explored, the FG and SG sub¬ 
systems have initially the same total mass and they both follow a 


Table 1 . Initial parameters for all models. 


Model id. 

N 

R h,FG 

R h,SG 

fb,0 

x g,0 

MPr5f03x3 

20,000 

5 

0.03 

3 

MPr5f03x5 

20,000 

5 

0.03 

5 

MPr5f03xl0 

20,000 

5 

0.03 

10 

MPr5f03x20 

20,000 

5 

0.03 

20 

MPrl0f03x3 

20,000 

10 

0.03 

3 

MPrl0f03x5 

20,000 

10 

0.03 

5 

MPrl0f03xl0 

20,000 

10 

0.03 

10 

MPrl0f03x20 

20,000 

10 

0.03 

20 

MPr5flx3 

20,000 

5 

0.1 

3 

MPr5flx5 

20,000 

5 

0.1 

5 

MPr5flxl0 

20,000 

5 

0.1 

10 

MPr5flx20 

20,000 

5 

0.1 

20 

MPr5flx3-20 

20,000 

5 

0.1 

3-20“ 

SPflx3 

20,000 

- 

0.1 

3 

SPflx5 

20,000 

- 

0.1 

5 

SPflxlO 

20,000 

- 

0.1 

10 

SPflx20 

20,000 

- 

0.1 

20 

SPflx3-20 

20,000 

- 

0.1 

3-20“ 

MPr5f03xl0n32k 

32,000 

5 

0.03 

10 

MPr5f03xl0n40k 

40,000 

5 

0.03 

10 

MPr5f03x 10n60k 

60,000 

5 

0.03 

10 


N = N s + No ; n is the total number of single, N s , and binary, ,\) nn , 
particles. 

fb 0 = Tv),i n / (N s + ATi n ) is the initial binary fraction. 

p h,FG is the ratio of the FG to the SG initial half-mass radii. 

K h,SG 

All the SP models refer to single-population systems. 
m gi 0 initial hardness parameter (see section 2 for definition). 

“ A uniform distribution in binding energy is assumed. 

King model (1966) density profile with central dimensionless po¬ 
tential Wo = 7; the SG subsystem, however, is initially concen¬ 
trated in the inner regions of the FG system as suggested by the re¬ 
sults ofD’Ercole et al. (2008). All the simulations start with equal- 
mass particles. We have explored two models here: one in which 
the initial ratio of the FG to SG half-mass radius (Rh,FG / Rb,SG) 
is equal to 5, and one with -Rh,FG/-Rh,SG = 10. In order to 
explore the differences between single-population and multiple- 
population clusters we have also run a number of simulations of 
single-population (SP) systems with the same number of particles 
and tidal radius as the multiple-population systems but with an in¬ 
ternal structure described by a single King model with Wo = 7. 

We include the effects of the host galaxy tidal field (modelled 
as a point mass) and assume that the cluster is initially tidally trun¬ 
cated; particles moving beyond a radius equal to two times the tidal 
radius are removed from the simulation. 

As far as the properties of the binary population are concerned 
we have explored the evolution of systems with values of the initial 
binary fraction, /b,o = Nbi n /(N S + A'bin) equal to 0.03 and 0.10 
(where N s is the number of single particles and A r bin is the number 
of binaries; we will refer to N = Noin + N s to indicate the sum of 
the total number of single and binary particles; the total number of 
particles in each simulation is equal to N to t = N s + 2IVbin). All 
the simulations start with N = 20, 000 (Ntot = 22, 000 for simu¬ 
lations with /b,o = 0.1, and Ntot = 20, 600 for simulations with 
fb,o = 0.03). We assume that FG and SG binaries have the same 
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Figure 1. Time evolution of the total number of binaries, A© n (t) (nor¬ 
malized to the total initial number of binaries, Af,j n (0)) for simulations 
MPr5fl and SPfl and different values of the hardness factor x g Q. Time is 
normalized to the initial half-mass relaxation time, T r h,o- 

initial binding energy distribution. In this paper we focus our at¬ 
tention on moderately hard binaries with initial global hardness pa¬ 
rameter, x g: o, ranging from 3 to 20; x gt o is defined as Eb/(mag P ) 
where Eb is the absolute value of the binary binding energy, and 
asp the 1-D velocity dispersion of all stars in the SP system. In or¬ 
der to better illustrate the dependence of the binary evolution on the 
hardness parameter, we have run a number of simulations each one 
including binaries with a single value of x g: o- We have also car¬ 
ried out one simulation including binaries with a uniform binding 
energy distribution between * gi o = 3 and x gl o = 20. 

The initial parameters of all models and the ids. used to 
identify them in the rest of this paper are summarized in Ta¬ 
ble |Tj Some of the simulations performed reach a deep core col¬ 
lapse phase; specifically the following systems reach the deep core 
collapse phase (MPr5f03x3, MPr5f03x5, MPrl0flx3, MPrl0flx5, 
MPrlOflxlO, MPrl0flx20, MPr5flx3, MPr5flx5, SPflx3) at 
t/ 7" r h,o equal to about, respectively, (8.4, 8.7, 4.9, 4.9, 5.7, 7.8, 7.3, 
8.3, 15.5). 

In order to provide an indication of the extent of the stochastic 
variations in the results obtained we have carried out 10 simula¬ 
tions with different random realizations of the initial conditions for 
the system MPr5f03x5 and 10 for the system MPr5f03x20. Three 
simulations with N = 32, 000, 40, 000 and 60, 000 and structural 
properties equal to those of MPr5f03xl0 have also been carried out 
to explore the dependence of the results on the number of particles. 

3 RESULTS 

3.1 Evolution of the global binary fraction 

Fig. □ shows the time evolution of the global binary frac¬ 
tion for the multiple-population systems MPr5flx3, MPr5flx5, 
MPr5flxl0, MPr5flx20, MPr5flx3-20, and the corresponding 
single-population systems SPflx3, SPflx5, SPflxlO, SPflx20, 
SPflx3-20. As shown by this figure, multiple-population systems 
are characterized by a general enhancement in the binary disruption 
compared to single-population clusters. As discussed in Vesperini 



Figure 2. Time evolution (with time normalized to the initial half-mass re¬ 
laxation time) of the total number of binaries (normalized to the total ini¬ 
tial number of binaries). At any given time the difference between 1 and 
the value shown by the solid line represents the fraction of the initial bi¬ 
nary population surviving in the cluster, the difference between the values 
shown by the solid line and the dashed line represents the fraction of the 
initial number of binaries ionized, and the dashed line shows the fraction 
of the initial number of binaries escaped from the cluster. Thick, red lines 
refer to the MPr5flx20 simulation; thin blue lines refer to the MPr5flx3 
simulation. 


et al. (2011), the differences in the structural properties between 
multiple-population and single-population clusters are responsible 
for the observed differences in the evolution of the global binary 
fraction. 

The role of different processes responsible for the decrease 
in the number of binaries is illustrated in Fig. [2] for the simula¬ 
tions MPr5flx3 and MPr5flx20. For the range of binding energies 
considered in this paper, binary ionization is the dominant process 
while binary escape play only a minor role. Ionization is more ef¬ 
ficient in the central high-density regions populated mainly by SG 
binaries and, as we will further discuss below, this leads to a sig¬ 
nificant preferential disruption of SG binaries. As expected, simu¬ 
lations with binaries with smaller values of x g ,o are characterized 
by a more rapid ionization of a larger fraction of binaries. 

Binary escape affects primarily FG binaries but the preferen¬ 
tial escape of FG binaries is far from being sufficient to balance 
the preferential ionization of SG binaries with the net result that 
the global fraction of SG binaries declines much more rapidly. Fig. 
[3] show the number of FG, SG and mixed binaries (hereafter we 
refer to binaries that, as a result of a component exchange event, 
have one SG and one FG component as mixed binaries) escaping 
for the MPr5flx3 and the MPr5flx20 simulations. This figure fur¬ 
ther illustrates the role of escape in determining the decrease of the 
global binary fraction and shows the extent to which SG, FG and 
mixed binaries are affected by this process. 

In Fig. [4] we show the time evolution of total number of SG 
and FG and mixed binaries. This figure clearly illustrates the pref¬ 
erential disruption of SG binaries in all the systems considered in 
this paper; it is also interesting to note that exchange events, partic¬ 
ularly for simulations with the harder binaries, lead to the formation 
of a non-negligible number of mixed binaries. 

A comparison of the panels of Fig. [4] showing the results for 
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Figure 3. Time evolution (time is normalized to the initial half-mass relax¬ 
ation time) of the total number of escaping binaries (normalized to the total 
initial number of binaries, A r |,j n (0)).Thick, red lines refer to the MPr5flx20 
simulation; thin blue lines refer to the MPr5flx3 simulation. Solid, dotted, 
and dashed lines refer, respectively, to FG, SG, and mixed binaries. 


the MPr5 and the MPrlO simulations, show that, as was to be ex¬ 
pected, systems starting with a more concentrated SG subsystem 
are characterized by a more efficient disruption of SG binaries. 
Even for the x g} o = 20 case, only about 6 % of the initial popula¬ 
tion of SG binaries survive after about 10r r h,o in the MPrl0f03x20 
simulation whereas in the less concentrated MPr5f03x20 system 
about 15 % of the initial SG binary population survives after the 
same number of initial half-mass relaxation times. Fig. [5] shows 
the time evolution of the ratio of the total number of SG to the 
total number of FG binaries and further illustrates the preferen¬ 
tial disruption of SG binaries. For the simulations MPr5f03x5 and 
MPr5f03x20, the plots in Fig.[4]and[5]show, respectively, the mean 
values of the number of binaries and the mean of the SG-to-FG 
binary number ratio along with the la variation about the mean as 
calculated from ten different realizations of each of these two initial 
conditions. 

In Fig. [6] we show the radial profile of the local hardness pa¬ 
rameter ®iocai(f), defined as the ratio of the binary binding energy 
to mcr(r) 2 where a(r) is the 1-D velocity dispersion measured at a 
distance r from the cluster centre. This parameter provides a more 
specific measure of binary hardness and of the expected dynamical 
fate of a binary as a function of the distance from the cluster centre. 
SG binaries are concentrated in the cluster inner regions where the 
larger velocity dispersion results in smaller values of xiocai(r) and 
a more rapid disruption. 

For binaries in the cluster outer regions (preferentially FG bi¬ 
naries), on the other hand, xiocal is large and binaries can easily 
survive. FG binaries initially orbiting in the cluster outer regions 
can be disrupted only later in the cluster evolution after migrating 
towards the cluster central regions as a result of the mass segrega¬ 
tion process. 

Fig[7] and [8] show the density of FG, SG, and mixed binary 
stars in the r-a;i OC ai(r) and the r-Eb planes for the simulation 
MPr5flx3-20. This figure illustrates the interplay between segrega¬ 
tion and binary evolution. As binaries segregate towards the inner 
regions they are disrupted as a result of encounters with single stars 


or other binaries. Harder binaries, however, can survive and evolve 
to become harder. SG binaries are those more affected by these 
evolutionary processes (ionization and hardening) as clearly shown 
also by the histograms of xi oca i(r) and Eb in Fig.[7|and[8] Fig. [7] 
and [8] also show that while segregation leads to a decrease in the 
number of binaries in the cluster outer regions (r > rh, where rh 
is the cluster half-mass radius), the distribution of binding energy 
of outer binaries is not evolving significantly. The effects of disrup¬ 
tion and hardening are instead clearly visible in the distribution of 
binding energy of inner ( r < rh) binaries. 

Finally, in Fig. [9] we show the time evolution of the to¬ 
tal number of FG and SG binaries for simulations with different 
initial number of particles (MPr5f03xl0n32k, MPr5f03xl0n40k, 
MPr5f03xl0n60k). Overall the results of our simulations show a 
small spread for systems with different number of particles. Since 
one of the mechanisms responsible for the decrease of the number 
of binaries is evaporation (although, as discussed above, it is not 
the dominant one) some spread is to be expected since the depen¬ 
dence of the evaporation time-scale on the number of particles dif¬ 
fers from that of the relaxation time (see e.g. Fukushige & Heggie 
2000, Baumgardt 2001). 

3.2 Radial variation of the binary fraction 

In this section we focus our attention on the radial profiles of the FG 
and SG binary fraction. As anticipated in the initial discussion in 
m the predictions concerning the different evolution of the global 
fraction of FG and SG binaries and the general preferential disrup¬ 
tion of SG binaries do not necessarily imply that the local binary 
fraction measured at different distances from the cluster centre fol¬ 
lows the same trend. 

In general, the details of the evolution of the radial variation 
of the SG-to-FG number ratio are determined by the complex inter¬ 
play among binary ionization, hardening, rate of segregation, and 
by differences in the initial spatial distribution of the FG and SG 
populations. 

While, as shown in the previous section, the overall fraction of 
SG binaries decreases more rapidly than that of the FG population, 
locally the binary population can be dominated by SG binaries. 

Fig. [TO] shows the time evolution of the radial profile of the 
binary fraction for all the binaries (upper panel), for FG binaries 
(middle panel) and for SG binaries (lower panel) for the simulation 
MPr5flx3-20. This figure illustrates the effect of segregation and 
binary disruption: segregation towards the cluster inner regions ini¬ 
tially leads to an enhancement of the central binary fraction, binary 
disruption offsets the effects of segregation and causes the binary 
fraction in the central regions to decrease. 

Fig. ED shows the cumulative (upper panels) and the differ¬ 
ential (lower panels) radial profiles of the SG-to-FG binary number 
ratio for the simulations MPr5flx5, MPr5flx20, and MPr5flx3-20. 
We point out that although the global FG binary fraction is larger 
than the global SG binary fraction (see Fig. 4), in the inner regions 
the binary population can still be dominated by SG binaries while 
in the outer regions, the binary population is dominated by FG bi¬ 
naries. 

Although our simulations are still focused on simplified mod¬ 
els and not aimed at a direct comparison with observational data, it 
is important to emphasize the difference between the global SG and 
FG binary fractions and the local values of these fractions measured 
at different distances from the cluster centre. Observational studies 
focusing on the binary fraction in a limited range of radial distances 
from the cluster centre might reveal a larger fraction of FG or SG 
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Figure 4. Time evolution of the total number of FG (solid lines), SG (dotted lines) and mixed (dashed lines) binaries for all the MPr5f 1 (first row), MPr5f03 
(second row), and MPrlOfOS (third row) simulations. Each column shows the results of simulations for a different value of :r g) o indicated at the top of the 
figure. For simulations MPr5f03x5 and MPr5f03x20 the lines show the mean of the number of binaries from the results of the ten simulations with different 
random realizations of the initial conditions and the shaded area show the =blcr. Time is normalized to the initial half-mass relaxation time. 




Figure 5. Time evolution (with time normalized to the initial half-mass 
relaxation time) of the SG-to-FG binary number ratio for MPr5f03 (upper 
left panel), MPrl0f03 (upper right panel), and MPr5fl (lower left panel) for 
different values of rr g) o (see lower right panel for colors and line styles cor¬ 
responding to the different values of :r g ,o)- For simulations MPr5f03x5 and 
MPr5f03x20 in the upper left panel, the lines show the mean of the SG-to- 
FG binary number ratio from the results of the 10 simulations with different 
random realizations of the initial conditions; bars showing the dzlcr are also 
shown at a few representative times. 


Figure 6. Radial profiles of initial local hardness parameter, £Ci OC alM (see 
§3.1 for definition) for the MPr5f03 (solid lines) and the MPrl0f03 (dashed 
lines) simulations. Different solid (and dashed) lines correspond (from top 
to bottom) to a: g) o = 20, 10, 5, 3. Radius is normalized to the half-mass 
radius, r h- 

binaries depending on the radial position of the region probed by 
observations. The interpretation of observational results in the con¬ 
text of the preferential disruption of SG binaries described in this 
paper must therefore take into account the radial variation of the 


© 0000 RAS, MNRAS 000, 000-000 




























6 J. Hong et al. 



-2 -1 0 1 -2 -1 0 1 -2 -1 0 1 0 0.05 0.1 0 0.05 0.1 0 0.05 0.1 

l°gio r / r h N bin /N bin (0) 


Figure 7. Time evolution of the density of SG (first column), FG (second column), and mixed (third column) binaries in the log(r/rh) — log(£i oca i(r)) 
plane (where x\ oca \(r) is the local hardness parameter; see §3.1 for definition) for the MPr5f 1x3-20 simulation. The last three columns of panels (on the right 
side of the figure) show the histograms of the distribution of log(xi oca i(r)) for SG (solid blue lines), FG (red lines and short vertical segments), and mixed 
binaries (dashed lines) for the whole system, for binaries within the half-mass radius and outside the half-mass radius. Data from three combined snapshots 
around the indicated time are used for these plots. 
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Figure 8. Time evolution of the density of SG (first column), FG (second column), and mixed (third column) binaries in the log(r/rh) — log(£^b/(-^b o)) 
plane (where E b is the binary binding energy and ( E b o) is the initial mean binding energy of all binaries for the MPr5f 1x3-20 simulation. The last three 
columns of panels (on the right-hand side of the figure) show the histograms of the distribution of log(£'b/(-^b o)) f° r SG (solid blue lines), FG (red lines 
and short vertical segments), and mixed binaries (dashed lines) for the whole system, for binaries within the half-mass radius and outside the half-mass radius. 
Data from three combined snapshots around the indicated time are used for these plots. 
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Figure 9. Time evolution of the total number of FG (thin lines) and SG 
(thick lines) binaries (normalized to the total initial number of binaries) for 
the model MPr5f03xl() and different values of the total number of particles, 
N (see legend in the figure). Time is normalized to the initial half-mass 
relaxation time. 



0.1 1.0 10.0 
r/r h 


Figure 10. Time evolution of the radial profile of the binary fraction for the 
model MPr5flx3-20 for all binaries (top panel), FG binaries (middle panel), 
and SG binaries (lower panel). Different lines in each panel correspond to 
the times indicated in the legend in the top panel. The radial profiles have 
been calculated with data from three combined snapshots around the indi¬ 
cated time. 


SG and FG binary fractions and the differences between local and 
global binary fraction of the two populations. 

While more realistic simulations are needed for a close com¬ 
parison with observational studies we have further addressed this 
issue here and explored the time evolution of the distance, f? g iob, 
from the cluster centre where the local FG-to-SG binary number 
ratio equals the global value of this ratio. Fig. [12] shows the time 
evolution of R g iob for all the simulations presented in this paper. 
.Rgiob does not vary significantly during the cluster evolution and 
it does not significantly depend on the SG initial concentration or 
the binary fraction. R g \ 0 h lies between approximately rh (for sim¬ 
ulations with high-binding energy binaries, a; g ,o > 5) and 2-2.5 rh 
(for simulations with low-binding energy binaries which are more 
efficiently disrupted at larger distances from the cluster centre). Ob¬ 
servations focusing on this range of radial distances should there¬ 
fore provide an approximate indication of the actual global value 
of the FG-to-SG binary number ratio. 


4 CONCLUSIONS 

In this paper we have presented the results of a survey of N- 
body simulations aimed at exploring the evolution of binary stars 
in multiple-population globular clusters. In this study we have ad¬ 
dressed a number of fundamental aspects of the dynamics of binary 
stars in multiple-population clusters; in subsequent studies we will 
further expand this investigation by including a spectrum of masses 
and explore how the presence of a spectrum of masses affects the 
interplay among mass segregation, binary ionization, hardening, 
and exchange interactions and the role of these processes in the 
dynamics of the binary population. 

Theoretical models for the formation of multiple-population 
clusters predict that SG stars form segregated in the cluster inner 
regions and a number of observational studies have shown that in 
several clusters SG stars are more centrally concentrated and still 
preserve some memory of the initial segregation predicted by the 
formation models. The results of the simulations presented in this 
paper show that the initial differences between the structural prop¬ 
erties of the SG and the FG populations may leave a fingerprint in 
the current properties of binary stars. 

We have shown that the SG binary disruption rate is larger 
than that of the FG binary population. This is a consequence of the 
more efficient binary disruption in the cluster central high-density 
regions populated mainly by SG binaries. In addition to disruption, 
escape also affects the number of binaries in a cluster and we have 
shown the relative role of these processes in driving the evolution 
of SG and FG binaries. 

The difference between the initial spatial distributions of the 
SG and the FG populations has also implications for the properties 
of the surviving binaries. Specifically, as a result of the more rapid 
hardening of surviving binaries in the cluster inner regions, the SG 
binding energy distribution is characterized by a larger fraction of 
binaries with large (more bound) binding energies. 

We have studied the evolution of the radial variation of the SG 
and FG binary fraction with the distance from the cluster centre. 
The radial variation of the binary fraction is the result of the initial 
differences between the spatial distribution of the SG and FG popu¬ 
lations and of the subsequent effects of binary segregation, disrup¬ 
tion, and escape. Although, as the cluster evolves, SG binaries are 
preferentially disrupted and the SG global binary fraction decreases 
more rapidly than that of the FG binary population, the local binary 
fraction in the cluster inner regions may be larger for SG binaries 
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Figure 11. Top panels: time evolution of the radial profile of the ratio of the cumulative radial distributions of SG to FG binaries, IVb sg(< r )/Nb FG (< r ) 
for the simulations MPr5flx5, MPr5flx20, and MPr5flx3-20. Bottom panels: time evolution of the radial profile of the SG-to-FG binary number ratio, 
Nb SG( r )/-^b fg ( r ) for the simulations MPr5flx5, MPr5flx20, and MPr5f 1x3-20; the horizontal lines indicate the global SG-to-FG binary number ratio . 
The radial profiles refer to the times shown in the legend in the lower panels. The radial profiles have been calculated with data from three combined snapshots 
around the indicated time (except for t/r T h o = 10 for which five snapshots were used). 
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Figure 12. Time evolution of Rg\ 0 b > the radius where the local SG-to-FG 
binary number ratio equals the global value of this ratio. Different symbols 
and colors correspond to different simulations as indicated in the legend. 
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while the outer regions are mostly populated by FG binaries. We 
have shown that the SG-to-FG binary number ratio measured at a 
distance from the cluster centre equal to about 1-2.5 r b (where rh 
is the cluster half-mass radius) is approximately equal to the global 


value of this ratio. Possible differences between SG and FG global 
and local binary fraction should be taken into account in the inter¬ 
pretation of observational studies aimed at exploring the properties 
of SG and FG binaries spanning limited range of distances from a 
cluster centre. 
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